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ABSTRACT 

We present a series of decaying turbulence simulations that represent a cluster-forming 
clump within a molecular cloud, investigating the role of magnetic fields on the forma- 
tion of potential star-forming cores. We present an exhaustive analysis of numerical 
data from these simulations that includes a compilation of all of the distributions of 
physical properties that characterize bound cores - including their masses, radii, mean 
densities, angular momenta, spins, magnetizations, and mass-to-flux ratios. We also 
present line maps of our models that can be compared with observations. Our simu- 
lations range between 5-30 Jeans masses of gas, and are representative of molecular 
cloud clumps with masses between 100 — IOQQMq. The field strengths in the bound 
cores that form tend to have the same ratio of gas pressure to magnetic pressure, /3, 
as the mean /? of the simulation. The cores have mass-to-flux ratios that are generally 
less than that of the original cloud, and so a cloud that is initially highly supercritical 
can produce cores that are slightly supercritical, similar to that seen by Zeeman mea- 
surements of molecular cloud cores. Clouds that are initially only slightly supercritical 
will instead collapse along the field lines into sheets, and the cores that form as these 
sheets fragment have a different distribution of masses than what is observed. The spin 
rates of these cores (wherein 20-40% of cores have fli// > 0.2) suggests that subse- 
quent fragmentation into multiple systems is likely. The sizes of the bound cores that 
are produced are typically 0.02-0.2 pc and have densities in the range 10^ — 10^ cm~'^ 
in agreement with observational surveys. Finally, our numerical data allow us to test 
theoretical models of the mass spectrum of cores, such as the turbulent fragmentation 
picture of Padoan-Nordlund. We find that while this model gets the shape of the core 
mass spectrum reasonably well, it fails to predict the peak mass in the core mass 
spectrum. 

Key words: MHD - turbulence - stars: formation - ISM: clouds - ISM: kinematics 
and dynamics 



1 INTRODUCTION 

Molecular clouds are highly dynamic environments and un- 
derstanding the effects that their internal, supersonic tur- 
bulent motions have on the process of gravitational col- 
lapse is critical towards understanding how stars form. The 
turbulent fragmentation picture of star formation posits 
that star formation is regulated by the turbulent motions 
of the molecular cloud; turbulent fluctuations can sup- 
port the cloud on large scale while also creating compres- 
sions that can become gravitationally unstable, wherein 
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protostellar cor es form after 1-2 free-fall times; 0.5-1 Myr 
(eg., reviews by l|Mac Low fc Klessenll2004l : lElmegreen et all 
'2000")). This is in stark contrast with the traditional 
paradigm of star formation, which supposes that a prestel- 
lar core in a magnetically subcritical cloud contracts qua- 
sistatically as the magnetic field slowly leaks out due to 
ambipolar diffus i on on time-scales of several Myr (e.g . 
'Mest el fc Spitze r'('l956'l:'Mouscho vias fc Paleologoul (Il98lf): 
iPaleologou fc Mouschoviaa (19811 ): IShu. Adams fc Lizanol 
l|l987h : [ciolek fc Basul (|200ll )). 



Our goal i n t his and the previous paper 



iTflley fc Pudrit j |2004 TP04 hereafter) is to under- 
stand the origin of star clusters, in particular the origin of 
the distribution of molecular cloud cores (and their physical 
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properties) out of which cluster stars are formed. Numerical 
investigations of turbulence within GMCs have shown 

that turbulent fragmentation is remarkably successful 

in produ cing clust ers of stars (IGammie fc Qstriker 1996 



Klessen e t al. 200^ 

iHeitsch et al] 



2001 



Padoan fc Nordlundl I2OO2I: iLi et a" 



199^ : 
200 ll: 



20051 : IVazguez-Semadeni et all bOOSbl ) 



.Qstriker et alj ' . . . 
20041 : iKlessen etla] 
Many of the ob- 
served statistical properties of cluster-forming regions in 
molecular clouds can be reproduced by these models, 
including: 

(i) the apparent relati on between internal v elocity 
and size (Burkert fc Bodenheimei 200C; 



2002 



Qstriker et al.l I2OO I: Ballcstcros -Paredes fc Mac Low! 
Qssenkopf fc Mac Low.. 2002 : .Hever fc Bruntll2004l ): 
(ii) the simil arity with th e obse rv ed distribution of 



core masses ('Klessen et al. 
(l2000l'): ..Balsara elal (2001 
KlessenI 



l2(301al lbl): 



, Klessen fc Burkert 

Reipurth fc Clarkd (I2OOII): iBate et al 
jallesteros-Paredes fc^fac Lowl ll2002l); 



Klessen & Burkert 
Heitsch et al. ( 2001) 
2001 



20031): lLietal.1 (l2003ll: 'Nordlund 
Qstr iker (2003) 



Bate et al] 



Padoan 
Clark & 



iGammie et al 

\2m4): iDelgado-Donate. Clarke fc Batel 
iDelgado-Donate et all (120041} : iKlessenI (|2003): 



('2003') 
Bonnc ll 
lf2004h':; 



Li et al 



(2004 [ l: iPadoan fc Nordlundl (l2004h : TP04: iBate fc Bonnel] 



1 20051): IJappsen et al.l ll2005h: iMartel et al 



iPadoan et al.l ( 2005l l: iBallesteros-Paredes et al.l 



- wherein the numerical data shares many of the key 
features of observed star-forming region, such as the 
(Salpeter) ~ —1.35 power-law relationship at large masses 
between number of cores of a given mass and the core 



mass and a turnover at low masses (iKramer et al.l 


19981: 


Motte et al.l 19981: Testi fc Sareend 19981: Ijiiina et all 


1999: 


Johnstone et al. 


IgOOO"; 'Luhman' '2000*; 'Luhman et al. 


2009: 


Johnstone et al. 


I2OOI: Briccno ct al. 2002: Luhman et al. 


20031: iReid fc Wilsonll2005l): 



(iii) the similarity of the co re mass function to the stellar ini- 
tial mass function (IMF) (dMiller fc Scald[l979l : IScaloll 19861 . 
1 19981 : iKroupallioO^ : IChabrie r"200?ri: 

(iv) the similarity of the a ngular momentum of cores wit h 



observed rotations o f cores dBurkert fc Boden 
iKlessen fc BurkertI lioool'): IGammie et al.l 



( 2004h : IJappsen fc Klessa] (|2004h : iLi et'al] l|2004h .TP04V 
which are consistent with obs ervations of v e locity - 
gradients in molecular clouds (iGoodman _et all I1993I : 
iBarranco fc Goodmanlll998l : ICaselli et al.ll2002al lbl) . and sev- 
eral orders of magnitude greater than the angular momenta 
of main- sequence stars; and 

(v) the moderate efficiency of star formation within dense 
clumps wherein a fraction of the gas in these models col- 
lects in th e cores, up to 30-40 % in models of decayiii g 
turb ulence ('Cammie et al' ("20031); IClark fc Bonnell' (2004); 
IVaza ucz-Scmadcni ot al. (2003):TP04 :lciark ot al. (200^)). 
Models with dr iven turbulence predict lower st ar formation 
efficiencies (e.g. lVazguez-Semadeni et al.l (|2003l )'). and it has 
been suggested that the star for mation needs to be spread 
out over s everal dynamical times ('Kr umholz fc McKeell2005l : 
iTan et al . 2006; Krumholz fc Tan 20^. 

Qur previous work on purely hydrodynamic turbu- 
lence in self-gravitating gas (TP04) found that only a 
few initial Jeans masses are needed to create a cluster of 
prestellar cores. This is because the shocks and compres- 



sions create local density enha ncements with much red uced 
local Jeans masses Jsee also iGoodwin etHI (|2004bl lalH): 
IClark fc Bormeiil (|2005l )l. The point here is that since the 
local Jeans mass is inversely proportional to the square root 
of the density, it can become significantly easier to form col- 
lapsing cores in these density enhancements. 

In the present paper, we extend our investigation to in- 
clude the role that magnetic fields can play in turbulent 
fragmentation. Magnetic fields can play a role in princi- 
ple because of the well known fact that they change the 
condition for gravitational stability. In a uniform magne- 
tized medium, the condition for collapse is most physically 
measured by the magnetic criticality parameter (or mass- 
to-fiux ratio) F = 2tv\/1JT,/B, where E is a surface density 
of the gas and B is its magnetic field strength. In situations 
in which the magnetic field is weak compared to gravity - 
the so-called magnetically supercritical case, P > 1 - grav- 
itational collapse is expected. The subcritical configuration 
(r < 1) has magnetic fields strong enough to support the 
cloud against gravitational collapse, and is thus expected 
to be stable. A condensation can collapse only if the mag- 
netic field can escape via ambipolar diffu sion, reducing the 



magne t ic fiux and thus increasing P (e.g . iMestel fc Spitzer 
l'l956t): lMouschovias fc Paleologoul (ligSll il. lLi fc Nakamura 



(2004) have demonstrated that it is possible for turbulence 
to enhance this diffusion rate, thus hastening the transition 
from subcritical to supercritical behaviour, but it is not clear 
if the cores formed from this type of process resemble GMC 
cores. 

Given that the unmagnetized simulations have had such 
success - what role does the magnetic field play? Based on 
the previous paragraph, it is obvious that in regions within 
the cloud wherein the field is too strong (and couples ideally 
to the gas), magnetic forces can prevent collapse from occur- 
ring. This has been c onfirmed in simulations of subcritica l 
MHD turbulence (e.g. iHeitsch. Mac Low fc KlessenI (|200lh : 



iQstriker. Stone fc Gammiel (|200ll ) ~ although it has been 



suggested that the magnetic field can leak out of conden- 
sations via turbulence-enhanced ambipolar diffusion in rea- 
sonable time-sca les, producing supercritic a l cores out of sub- 
critical regions (iFatuzzo fc Adams! I2OO2I: I Kim fc Diamond! 
I2OO2! : IZweibell I2OO2I : iLi fc Nakamura! |2004| ). If the initial 



cloud is magnetically supercritical, collapse will occur. An- 
other important consideration that is not well handled by 
idealized analytical models of clouds is the fact that ge- 
ometry of the magnetic field is important since the field 
does not resist motions parallel to it. A third factor that 
makes magnetic fields potentially interesting is that any ro- 
tation within the cores can generate torsional Alfven waves 
that can carry away excess angular momentum and possibly 



that will ultimately form ("Mouschovias fc Paleologod 


197g, 


1980; 


Blandford & Payne 


1982; Pudritz fc Norman] 
fc Mouschoviadll994l). Do 


1983; 


Lovelace et al.j|l99ll: iBasu 


these 



important effects influence the efficiency of forming bound 
cores in turbulence? 

In this paper, we examine the effect of a magnetic field 
in addition to turbulence in self-gravitating, cluster-forming 
clumps within GMCs. We perform an exhaustive examina- 
tion of the contributions to the virial equation from the ther- 
mal, kinetic, gravitational, magnetic and surface terms of the 
fiuctuations that are produced by the turbulence. We use the 
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virial equation to identify the objects that are bound or col- 
lapsing, which we identify as 'cores'. The cores are confined 
through a combination of thermal, dynamic and magnetic 
surface pressure, and gravity. The virial equation provides a 
useful guide to determine the forces acting on the cores. In a 
turbulent environment, the surface press ure terms can have 
a sign i ficant role in confining th e core (|McKee fc Zweibel 
jim'i ): [Sallesteros-Paredes et all ([l99?),TP04). An appar- 
ent equipartition, at least in a statistical sense, develops 
between the internal kineti c ener gy and the gravitational 
energy (TP04 lKlessen et all l|2005h ). 

We find that strong initial magnetic fields (F ~ 1 — 5) 
can indeed have a large effect on the resulting evolution 
of the cloud and result in collapse to large sheets. There 
is little evidence for such prevalent structures from obser- 
vations however, and we use this to constrain the physical 
parameters of our simulations to uncover the conditions that 
are most favourable to the formation of observed star clus- 
ters. The simulations that initially have large F appear to 
do much better at reproducing the observed trends seen in 
molecular cloud mass functions. These simulations can still 
produce cores with local values of F that are close to crit- 
ical, as the cores are produced when the turbulence breaks 
up the fluid into parcels with smaller masses and magnetic 
fluxes than the original computational domain. We also flnd 
good agreement between our core properties with ammo- 
nia surveys of core s in clustered star formation regions (e.g. 
iJiiina et al.l |l99i)). 

We also use our numerical data to compute the an- 
gular momenta and spin rates for the cores in our simu- 
lations. We flnd that the angular momenta arc on the order 
of io^^-^3 Q-Q^g-i^ compara ble to measurements of GMC 



cores IIGoodman et al. I I1993I : [Barranco fc Goodman! 1199^ : 
iGaselh et al.ll2002al lbh. The distribution of core spins that 
we find suggests to us that the binary frequency of stars 
that form within these objects, must be high. 

This work, wherein we compare our results with a wide 
variety of observations of cloud cores, shows that the strong 
magnetization of cores that is often observed is the result of 
local compression by the turbulence, and is not characteris- 
tic of the entire volume of the cloud. This and other results 
suggest that it is turbulence, and not the wide-spread infiu- 
ence of a powerful cloud magnetic field, that is central to the 
origin of core formation and core properites and ultimately, 
the origin of the IMF. 

The setup of our simulations is described in Section 2, 
and focuses mainly on the differences from the simulations in 
TP04. In Section 3, we describe the overall results of the sim- 
ulations, noting especially the morphological structure that 
is produced. In Section 4, we examine the dynamical state 
of the cores and condensations that form. We present, in 
Section 5, the distribution functions of many physical prop- 
erties of cores, including their masses, angular momenta, 
spins, radii, average densities, magnetizations and mass-to- 
flux ratios. We compare these with the observations and use 
them to test theoretical models of turbulent fragmentation. 




10 
nj 



100 



Figure 1. Initial conditions for the simulations presented in this 
paper. The various data points correspond to the values in Table 
[1] Superimposed are lines of constant magnetic critical number 

r. 



2 SIMULATIONS 
2.1 Initial Conditions 

The simulations presented here are set up in a very similar 
manner to those in TP04, and a re described in more det ail 
there. We use the zeus-mp code (|Stone fc Normanlll992al lbl: 
iNormanI l2000l l. a magnetohydrodynamic fluid dynamics 
code made available by NCSA. MHD f orces are calcu- 
lated using the Method of Characteristics (|Hawlev fc Stonj 
I1995D : self-gravity of the fluid is calculated using the FFTW 
fast fourier transform libraries (|Frigo fc Johnsonlll998l ). Our 
boundary conditions are set to a periodic grid, and represent 
the centre of a GMC. The calculations presented here were 
performed on a Compaq AlphaServer SC40 at the SHAR- 
CNET McMaster University site. These simulations have a 
resolution of 256^ grid cells, and were typically simulated 
using 8 processors. As the cores we wish to simulate have 
densities of 10'^ — 10^ cm~'^, an isothermal equation of state 
is appropriate ()Havashilll96g : lLarsonlll969l l. We have grav- 
ity active from the beginning of our simulations, but this 
does not have a signiflcant effect until the turbulence has 
significantly decayed. 

We start our simulations with an initially uniform den- 
sity field. We initiate the turbulence by creating a set of 
plane waves in Fourier space, with relative amplitudes drawn 
from either a 7 — —5/3 ID Kolmogorov spectrum (marked 
'K' in Table [ij or a 7 = —2 ID Burgers spectrum (marked 
'B' in Table [l| such that the velocity fiuctuations as a func- 
tion of 3D wavenumber k = 2tt/X scale as oc k"''^ . We 
remove waves longer than 1/8 the length of the box in order 
to minimize the effects of the periodicity of the box on the 
kinematics. We also imposed an exponential cutoff on short 
wavelengths, to eliminate fluctuations in the initial condi- 
tions on scales we are not able to properly resolve. This 
spectrum is illustrated in a later flgure. The initial veloc- 
ity field is purely solenoidal, but otherwise each plane wave 
is given a random phase and direction. We do not provide 
any forcing of the turbulence at later timesteps; the initial 
kinetic energy is allowed to decay freely. After ~ 1 shock- 
crossing time, the turbulence is fully developed. 
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Run 


Spectrum 


rij 




L/pc 


M\ 


/3 


MaX 


r 


B5b 


K 


4.6 


105.1 


0.32 


5.0 


0.9 


4.7 


4.9 


B5c 


K 


4.6 


105.1 


0.32 


5.0 


5.0 


11.2 


11.6 


B5d 


K 


4.6 


105.1 


0.32 


5.0 


10.0 


15.8 


16.4 


C5c 


K 


7.5 


453.9 


1.0 


5.0 


5.0 


11.2 


13.7 


C5d 


K 


7.5 


453.9 


1.0 


5.0 


10.0 


15.8 


19.3 


C5e 


K 


7.5 


453.9 


1.0 


5.0 


50.0 


35.4 


43.2 


D5a 


K 


12.0 


623.7 


1.0 


5.0 


0.1 


1.6 


2.3 


D5b 


K 


12.0 


623.7 


1.0 


5.0 


1.0 


5.0 


7.1 


D5c 


K 


12.0 


623.7 


1.0 


5.0 


3.9 


9.9 


14.1 


E14b 


B 


27.5 


1086.3 


1.0 


14.1 


1.0 


14.1 


9.4 



t M = t;RMs/cs 
t Ma = v-RMs/vA 



Table 1. Initial conditions for the simulations presented in this paper. The first letter specifies the number of Jeans masses - 'B' for 
nj = 4.6, 'C for 7.5, 'D' for 12.0, and 'E' for 27.5. The number represents the initial RMS thermal Mach number of the simulation, 
'jVf'. The final letter denotes the mean /3 of the simulation, where f} is the ratio of thermal pressure to magnetic pressure, 'a' represents 
P = 0.1, 'b' represents /3 ft: 1, 'c' represents /3 4 — 5, 'd' represents /3 = 10, and 'e' represents f) = 50. 



The initial magnetic field is uniform in strength and 
direction. We characterize its initial strength through the 
ratio of thermal energy to magnetic energy, 

/3 = AT,pcl/B^ (1) 

We select a rang e of j3 be t ween 0.1 and 50.0 (the Zeeman 
measurements of ICrutcherl (|l999l ) suggest f5 « 0.1); the ini- 
tial parameter space of {nj,P) is plotted in Fig. ([T|. 

The list of our various model parameters is given in 
Table [1] where models with initial Jeans number n,j — 
4.6, 7.5, 12.0 and 27.5 are designated with letters 'B','C','D' 
or 'E' respectively. The following number is the RMS Mach 
number of the initial turbulence spectrum; all of our simula- 
tions save one (E14b) were run with a turbulent Mach num- 
ber M = WRMs/ca = 5. We chose this value because experi- 
ence with our hydrodynamic studies (which featured models 
similar to our 'B' models) showed this to be a good regime 
for vigorous turbulent fragmentation into a well-populated 
CMF. The final, lower-case letter gives an indication of the 
strength of the initial magnetic field, with /3 ~ 0.1, 1, 5, 10, 50 
indicated respectively by a,b,c,d,e. 

The turbulent amplitude can also be viewed in terms 
of the Alfven Mach number. Ma = vrms/va where va ~ 
B/^/lAnp) = Cs/Vl^ is the signal speed of MHD transverse 
waves. The values of this number are given by 



density p, size L and magnetic field strength B - can be 
rewritten as 



Ma = x/STW 



(2) 



which for M = 5 and our range of (3 gives 1.6 < Ma < 35. 
Thus, all the simulations in this paper were performed in the 
mildly to strongly super-Alfvenic regime. The initial level of 
turbulent kinetic energy is greater than the gravitational 
self-energy of the cloud, but as the turbulence decays the 
gravitational force will begin to dominate the dynamics. 

Another important characteristic of the magnetic field 
that is useful in determining the behaviour of the fiuid is 
whether it is magnetically supercritical or subcritical. As 
already noted, a magnetically supercritical fluid has suffi- 
cient mass for gravity to overwhelm magnetic support and 
go into collapse; a magnetically subcritical fluid does not 
have enough mass, and thus will not collapse gravitation- 
ally. The critical mass to flux ratio - which is the ratio of 
the gravitational to magnetic energies in a parcel of gas of 



r = 



2-!vG^^^pL 



B 



3.12/3 



1/2 1/3 



(3) 



(IBasu fc Ciolekll2004l ). We plot our initial simulation mod- 
els in our 2-D parameter space (/3,nj) in Fig. ((T|, where 
the solid lines correspond to constant values of F. The value 
F — 1.0 is the dividing line between supercritical and sub- 
critical magnetic fields and collapse will occur for F > 1. 
The scaling of the critical parameter in Equation ((3} with 
our two fundamental parameters, /3 and nj, shows that 
strongly magnetized cores (/3 <^ 1) can still be supercritical 
if the initial number of Jeans masses nj is sufficiently large 
(n.7 > 0.36/3"^''^). Thus, sufficient gravity (i.e. large nj) can 
still leave a strongly magnetized cloud supercritical. 

Our simulations were run until we could not properly 
resolve the collapse everywhere on the grid, a process that 
generally took about 3 fiow-crossing times to occur. When 
scaled to the initial values of L in Table [T] this is ~0.5 Myr, 
consistent with observati ons that star formation generally 
lasts less than a few Myr l|Hartmannll200ll : iHartmann et al.l 
I2OOII ) . We used the local Jeans length as our criterion for this 
resolution, calculated for each cell using the density of that 
cell. The simulation stopped when the local Jeans length 
of a cell somewhere on the grid was less than four zones 
jTruelove et al.lll997l ) (this is equivalent to an increase in 
density of 4148. 6nJ^^^ over the initial density of the simu- 
lation; for our simulation values of nj G [4.6, 27.5] this cor- 
responds to final-to-initial density ratios of 450-1500). This 
criterion was established for simulations in the absence of 
magnetic fields; Hcitsch ct al. [2001]) have suggested a more 
stringent criterion of resolving the Jeans length by at least 
six zones when magnetic fields are included. However, our re- 
sults should not be contaminated by artiflcial fragmentation 
due to our less stringent criterion as we stop the simulations 
immediately upon the violation of this criterion, and the col- 
lapse of the region that violates this condition is happening 
on much faster time-scales than the dynamics anywhere else 
on the computational grid. As a result, any instabilities that 
might result due to violation of the more stringent Heitsch 
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et al. criterion will not have a chance to propagate to other 
cells before the simulation ends. 



2.2 Comparison of Units with Observations 

These simulations are scale-free, and can thus, in theory at 
least, be scaled to any size or mass. However, it is not nec- 
essarily physical to do so. As in TP04 we chose the initial 
condi tions from among the range of clumps and cores seen 
in the lLada et al.l (|l99ll ) survey; these are cores with masses 
between 10 — SOOM© and sizes between 0.05 — 0.5 pc; they 
contain 0.5 — 30 Jeans masses, assuming they are at a tem- 
perature of 20 K. This is roughly the range we expect our 
simulations to scale while still being representative of clumps 
and cores; at larger scales, we would expect clouds to have 
many more Jeans masses. 

The mass, size and magnetic field of our simulations 
scale as 



B 



119n 



2/3 



16.8nj/3"^/^ 



(- 



T 
20 K 



mtot 



V lOOmo 



(4) 
(5) 



2.3 Identification of Cores 



We identify individual condensations in our simula- 
tions through the application of a watershed a lgorit hm 
l|Vincent fc Soillelll99ll : iMangan fc Whit akeij 1 19991 ) that we 
developed in TP04. This algorithm efflciently identifies cores 
by identifying the local gradient vector at each point. A path 
is then traced, from zone to zone, following the local gradi- 
ent vector of each cell until a local maximum is reached. All 
zones along this path are assigned to a core marked by that 
local maximum. See TP04 for detail s, including a compar i- 
son with the CLUMPFIND algorithm of lWiUiams et al.l (|l994h . 

As the algorithm by default is sensitive to small fluctua- 
tions, it tends to break up larger cores into multiple smaller 
cores. We reduce the influence of this effect on our results 
by applying the core-flnding algorithm to the density field, 
smoothed with a Gaussian kernel with a width of 3 zones. 
This significantly reduces the number of small-scale cores 
that we find. 



3 CLOUD STRUCTURE 

3.1 Sheets, Filaments and Cores 

The structures that form in our simulations are shown in 
Figs. I2I5I and are delineated by isodensity contours and 
magnetic field lines. We find that there is a clear trend 
in the results for the fluid to flow preferentially along 
the m agnetic field lines, rather than perpendicu lar to the 
fields l|Padoan fc Nordlundl [l999l : IlI et all I2004D . This ef- 
fect is naturally more significant in the simulation runs 
with stronger magnetic fields (i.e. lower /3 and, in partic- 
ular, low r). The simulations with a mean magnetic field 
strength that is closer to the critical magnetic field strength 
show that fluid motions perpendicular to the mean mag- 
netic field are reduced, while motions parallel to the mean 
field are not significantly affected. In this regime the ma- 
terial forms a sheet that then fragments. The formation 



of a sheet is a well-known behaviour of subcri tical cores 
(|Mouschovias fc Spitzerlll976l : lHeitsch et aLlbOOll ). although 
such clouds would not be able to collaps e in the absence 
of ambipolar diffusion ( , Balsara et al. 1200 ll : iLi fc Nakamural 



I2OO4I : IVazguez-Semadeni et al. 2005b ). However, we have a 
sheet forming in a marginally supercritical cloud. We still 
show in Sections and ((5)l that the cores that form within 
this sheet have different statistical properties from observed 
cores and are unlikely to be representative of star-forming 
regions. 

The fluid in the low /3 simulations initially clumps to- 
gether locally, forming many small condensations. These 
condensations are gravitationally attracted to each other, 
but the component of the motion perpendicular to the field 
lines is not generally strong enough to exert a significant 
perpendicular velocity unless the magnetic field is strongly 
subcritical. These fluid condensations will instead flow along 
the field lines, until the gravitational force parallel to the 
field is reduced as a sheet forms. As this sheet builds up in 
mass, the gravitational forces transverse to the field become 
strong enough to cause local collapse in several parts of the 
sheet, that then experience the typical runaway growth of 
a gravitationally collapsing object. This results in several 
massive cores at the point where we are forced to stop the 
simulation due to violation of the Jeans criterion, as opposed 
to the non-magnetic and strongly subcritical runs where one 
object tends to runaway before any other condensations are 
significantly evolved. 

Furthermore, the separation of the simulation into a 
high-density sheet and low-density envelope divides the con- 
densations found into two distinct categories: those with a 
high average density (in the sheet), and those with a low 
average density (in the atmosphere). We show in Section 
((5| that these two sets have distinct properties; as a result, 
these sheets do not look like places where clusters of stars 
will form. 

3.2 Magnetic Field Structure in Turbulent Clouds 

The magnetic field in a turbulent molecular cloud can be- 
come quite complex, as shown for two of our runs, C5d 
and E14b, in Figs. [3] and [S] Oblique shocks generate vor- 
tical motions in the gas which in turn drive torsional Alfven 
waves that propagate through the fluid, resulting in the 
twisting and tangled structures we see in these two fig- 
ures. These propagating torsional waves create regions of 
helical magnetic field as described in theoretical work by 
iFiege fc Pudrit3 l|2000l ) and which may have be en observed 
in Zeeman observations (|Robishaw fc Heiles! '2005l . The field 
lines can get quite twisted and tangled, even with a relatively 
strong magnetic field, as shown in Figs. I2I6I The net direc- 
tion of the initial field is readily apparent, as magnetic flux 
is conserved to numerical accuracy. 

We also show a series of close-up views around the third 
most massive core (zooming in by factors of 2 and 8/3) in 
Run E14b in Fig.|6l which is typical of the mid-mass isolated 
bound cores produced in these simulations. We chose this 
particular core to show the details around cores that have 
not begun to experience a runaway collapse. We can also see 
the proximity of this core to two smaller, nearby cores. The 
field lines around the target core are not strongly tangled. 
As there is less turbulent energy on small scales (by virtue of 
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Figure 2. Isodensity surfaces and magnetic field configuration for the entire volume of Run B5c. The edge of the grid has a length of 
0.32 pc when scaled to our default values in Table [T] 



the Kolmogorov spectrum), it is perhaps not too surprising 
that the fluid motions are unable to tangle up the field to a 
significant degree in these smaller cores. No pinching of the 
field due to gravitational collapse is seen, as the mass-to-fiux 
ratio for this core is F = 0.9. Examining the virial terms we 
can see that this core is bound primarily due to the presence 
of a surface pressure; gravity is only responsible for 25% of 
the virial confinement. This core contains ~ 40000 zones, so 
it is well-resolved; it has a density ratio between centre and 
surface of ~ 9, with a central density (when scaled to the 
initial simulation values in Table [TJ corresponding to 8 x 
10^ cm~^. We note that this is not the core that underwent 
runaway collapse that stopped the simulation, and we do see 
evidence of focusing in the magnetic field lines there, but it 
is remarkable that these cores have grown so large without 
significant motions transverse to the initial field. 



Any rotation of cores like these would be expected to 
generate torsional waves that extract their angular momen- 
tum. We see some evidence for this in the twisting of the 
field lines beneath the core even though this core has not 
had enough time to fully rotate. Moreover, in this particu- 
lar core, we find that the angular momentum vector L and 
mean magnetic field B are misaligned by 45 degrees. We 
examine this in greater detail in a subsequent paper, but 
note that similar studies have suggested that cores do not 
strongly align with either loc al or global magnetic field di- 
rections l|Gammie et al.|[2003t ). 



3.3 Line Profiles 

We present simulated line profiles along with a column den- 
sity contour map for Run E14b in Fig. [T] The top row uses 
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a line-of-sight parallel to the mean field. The bottom row 
uses a line-of-sight perpendicular to the mean field, with 
the mean field in the horizontal direction. These line profiles 
(shown in the right-hand panels of Fig. (Tj are calculated in 
the same way as TP04, by summing along each line-of-sight 
the density value of each zone multiplied by a Gaussian ex- 
ponential term centered on the line-of-sight velocity of that 
zone; the width of the Gaussian is determined by the amount 
of thermal broadening expected for a gas at a temperature 
of 20K. These profiles represent optically thin emission lines 
that emit over the entire range of densities that we have in 
our simulations, typically 10^ — 10^ cm~"^ (ammonia or CS, 
for instance) . The amplitudes of the lines are scaled relative 
to the highest-amphtude line. The abscissa of the line plots 
are in units of the sound speed. 



orders of magnitude, but the range in column density in the 
projection along the mean field is only a factor of 14 due 
to some preference for collapse along the field lines. As a 
result, the line profiles all tend to have similar amplitudes. 
The profiles become noticeably more asymmetric near the 
regions of highest column density as material streams on to 
the large cores through channeling by the magnetic field. 
This is not seen in the projection perpendicular to the mean 
field, which has a much wider variation in column density 
(~ 60) and less motion transverse to the magnetic field. As a 
result, the lin e profiles have a strong Gaus si an character, as 
observe d (e.g. Falgarone fc Phillip g (Il990fl : Schneider et al.l 
(199d ); lGoodman et al.l (|1998 ): .Park et all |20o1 VK 



The range of densities in Run E14b varies over four 




Run 




/3 


r 


SFE 


B5b 


4.6 


0.9 


4.9 


0.54 


B5c 


4.6 


5.0 


11.6 


0.76 


B5d 


4.6 


10.0 


16.4 


0.76 


C5c 


7.5 


5.0 


13.7 


0.64 


C5d 


7.5 


10.0 


19.3 


0.64 


C5e 


7.5 


50.0 


43.2 


0.46 


D5a 


12.0 


0.1 


2.3 


0.52 


D5b 


12.0 


1.0 


7.1 


0.56 


D5c 


12.0 


3.9 


14.1 


0.47 


E14b 


27.5 


1.0 


9.4 


0.14 



Table 2. The star formation efficiencies for each of the simula- 
tions at the time the Jeans condition is violated. This measures 
the fraction of the total mass that is contained within gravita- 
tionally bound cores. 



4 FORMATION AND EVOLUTION OF CORES 



4.1 Energy Balance 

We found in TP04 that an excellent understanding of the 
dynamical state of the cores was provided by the indi- 
vidual terms in the virial equation. We apply the virial 
equation using Eulerian coordinat es to the case of a mag- 
netized, self-gravitating fluid as ijMcKee fc Zweibel Il992l : 
iBallesteros-Paredes et al.lll999l) : 



ll+l [ PJ-'v -dS ^U + K + W + S + M + F=l-I' (6) 
2 A J g Z 
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Figure 5. Isodensity surfaces and magnetic field configuration for the entire volume of Run E14b. The edge of the grid has a length of 
1.0 pc when scaled to our default values in Table [T] 



where 



I = 


/ pr^dV 
Jv 


(7) 


u = 


3 / PdV 
Jv 


(8) 


K = 


[ pv'^dV 
Jv 


(9) 


W = 


- [ pr- V^dV 
Jv 


(10) 


s = 


- / [Pr + r- (pvv)] ■ dS 
Js 


(11) 


M = 


i- f B'dV 
47r Jv 


(12) 


F = 


f r-TM-dS 

Js 


(13) 



where v is the velocity, p is the density, P = pCg is the 
isothermal pressure, $ the gravitational potential, B the 
magnetic field, and T the Maxwell stress-energy tensor. The 
terms /, U, K, W, S, M, F indicate respectively the moment 
of inertia of the fluctuation, the internal thermal energy, the 
internal kinetic energy, the gravitational energy, the sum of 
the thermal surface pressure and dynamic surface pressure, 
the internal magnetic energy, and the magnetic surface pres- 
sure. 

The volume integrals are calculated over all cells iden- 
tified as part of the fluctuation. The surface integrals are 
calculated from the fluxes through every zone identifled as 
being on the boundary between one fluctuation and the next. 

Of these terms, the internal thermal, kinetic and mag- 
netic energies will always be positive; the other terms can 
be either positive or negative. We plot the sum of the sur- 
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Figure 6. Zooming in on the tliird most massive core in Run E14b. The top-left panel shows the entire simulation, the top-right panel 
zooms in by a factor of 2, the lower-left panel zooms in by a further factor of 8/3, and the lower-right panel has the same magnification 
as the lower-left, but rotated by about 90 degrees to show the field line structure. 



face terms S and F against the gravitational term W, each 
normalized to the sum of the internal terms, in Fig. [8] As 
in TP04, we label the fluctuations for which /' < "bound 
cores" , and fluctuations for which 7' > "unbound fluctua- 
tions" . 

To interpret Fig. [S] we recall that the virial equation 
describes a steady state when / = in Eq. [G] This corre- 
sponds to the straight line in Fig. |5] A core that falls on 
this line is in virial balance, with the support provided by 
the internal kinetic energy K, thermal energy U and mag- 
netic energy M balancing the conflnement provided by the 
gravitational energy W, the magnetic pressure and tension 
on the surface F, and the surface thermal-fram pressure S. 
A core will collapse if it lies below the line in Fig. [S] 

The highly supercritical simulations such as Run C5e 



show a virial plot that is very similar to the plots we found 
in TP04, where there were typically one or a few cores 
that were strongly gravitationally dominant, and the bulk 
of the cores either lying close to the line 7' = or with 
W/ (7C 4- [/) ~ (in TP04, there was no magnetic field, and 
thus M=0). Most of the rest of the bound cores are only 
marginally virial unstable, while the unbound cores can be 
quite far from virial stability due to significant internal pres- 
sure overwhelming the confining forces. 

We see in Fig. |S] that the surface terms have a more 
significant contribution to the virial equation as the mean 
magnetic fleld becomes stronger. In particular, the two runs 
with the smallest F, Run B5b and Run D5a, have surface 
terms that act to confine some of the cores that are twice as 
strong as the internal kinetic, thermal and magnetic energies 
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Figure 7. Simulated column density contours and line maps of Run E14b, the simulation presented in this paper that most resembles 
a star-forming cluster. The top row contains a projection along the mean field axis (left panel) and the accompanying line map (right 
panel). The bottom row contains a projection perpendicular to the mean field axis, with the mean field in the horizontal direction. See 
Fig. [5] for a 3D projection of this simulation. The abscissas of the individual line map plots are in units of the thermal sound speed. 



that act to provide support. The formation of a sheet in these 
simulations leads to an enhanced thermal surface pressure 
(pCg) on the cores within it as the density in the sheet is 
higher, and the motions of the fluid along the field lines 
leads to an enhanced dynamic pressure out of the sheet. 

The simulations with intermediate values of the crit- 
ical parameter F have many bound cores that are signifi- 
cantly out of equilibrium, compared to the highly supercrit- 
ical runs. These cores are strongly confined both by surface 
effects and by gravity. 



4.2 Kinetic energy evolution 

We plot the kinetic energy Et (normalized to their ini- 
tial values ai t — 0, Eko) of the simulations as a func- 
tion of time (in units of the fiow-crossing time tflow, 
the time it takes for a shock moving at the RMS 
velocity to cross the simulation domain) in Fig. (|9]). 
Our results agree with previously reported results that 
the kinetic energy deca ys as , marked b y the solid 



the kinetic energy deca ys as f , marked b y tne solid 
straight hue in Fig. (l9t IIMac Low et al.|[l998l: Istone et all 



19981: iBis kamp fc Mii Ue^ Il999l: IP adoan fc Nordlundl llQQgl - 
Ostriker et al. 2001; C ho fc Laza rian 2003). In all cases, col- 
lapse occurs after ^ 1.5 — 2 tflow, accompanied by an increase 



in the kinetic energy due to the presence of supersonic infall 
to the most massive cores. At this point the gravitational en- 
ergy of the simulation dominates over the kinetic, thermal 
and magnetic energies. 

The energy spectrum for Run E14b is also plotted in 
Fig. (|9]) at the time of the analysis. This run began with 
a power spectrum, but the kinetic energy evolved to 
a Kolmogorov spectrum over the inertial range, although 
there is an excess of energy on large scales. In comparison 
with the initial state, we see that substantial damping from 
numerical dissipation in the smallest scale (largest k) modes 
has taken place by the end of the simulation. 



4.3 Effect of fields on star formation efflciency 

In TP04, we estimated an upper limit to the star formation 
efficiency (SFE) of the simulations, where the SFE is defined 
as the ratio of the mass contained within bound cores to the 
total mass in the simulation. The SFE was 40-50 per cent 
for simulations with 4.6 Jeans masses, and higher for simu- 
lations with fe wer Jeans masses . This is comparable to the 
estimates that IClark fc Bomieiil (|2004D found for unmagne- 
tized, unbound clouds. 

The SFE for each of the simulations discussed in this 
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Figure 8. The relationship between the various terms in the virial equation in each of the runs shown in Table [T] The models are 
arranged in order of decreasing V. The surface terms are plotted against the gravitational term, each normalized to the sum of the 
internal thermal, kinetic and magnetic energy. The solid line in each plot is /' = 0. 



paper are presented in Table [2] IVazguez-Semadeni et al] 
l|2005al ) noted that the SFE systematically decreased as the 
magnetic field strength increased. We see some evidence for 
higher SFE with larger F, supporting this trend, but with 
notable outliers for Run E14b (with notably small SFE, pos- 
sibly related to the greate r levels of init i al tur bulence in this 
simulation) and Run C5e. IPassot et alj (Il995l ) report a sim- 
ilar trend for higher SFE at intermediate field strengths, 
which they attribute to magnetic braking acting to reduce 
the local shear that could otherwise disrupt a core. The mag- 



netic field channels the flow into the sheet; when the sheet 
fragments gravitationally, a much larger fraction of the mass 
is accumulated in self-gravitating cores. Only at very large 
field strengths does the field retard the collapse in these 
sheets. 

Our definition for the SFE implicitly assumes that a sig- 
nificant fraction of the gas in a core will collapse as a part 
of the protostar. However, there are several physical pro- 
cesses that can occur in these later stages that our model 
does not take into account, some or all of which will oc- 
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Figure 8 - continued 
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Figure 9. (Left) Kinetic energy as a function of time for tlie different runs presented in this paper. (Right) The energy spectrum of Run 
E14b, multiplied by The dashed curve represents the initial spectrum of the turbulence and the solid line the resulting spectrum 

at the time the Truelove conditions was violated, 7 flow-crossing times later. 



cur in a cluster-forming region like that represented by our 
simulations. Ultraviolet radiation and ionization from the 
most massive protostars will heat the gas, increasing the 
Jeans mass and thermal support and reducing the accretion 
rate (jPranco et al.|[l993 ). Stellar winds can blow accretin g 
material away from the protostar l| Adams fc Fatuzz"olll996l l. 



Magnetically-driven outflows can return gas to the cloud 
(iPudritz fc Norman|[l98^ : IShu et al.|[l98i : iKrumholz et alj 
l2005h . All of these processes will serve to decrease the SFE 
we estimate. However, most of the material that was orig- 
inally in the core probably ends up in a star via accretion 
through a disc. 
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Run 



(measured) Mj*" (predicted) 



B5b 


4.6 


4.9 


6.3 X 10"* 




0.13 


B5c 


4.6 


11.6 


7.9 X 10-4 




0.020 


B5d 


4.6 


16.4 


1.3 X 10-3 


7.9 


X 10-3 


C5c 


7.5 


13.7 


2.5 X 10-3 




0.020 


C5d 


7.5 


19.3 


1.6 X 10-3 


7.9 


X 10-3 


C5e 


7.5 


43.2 


4.0 X 10-3 


1.0 


X 10-3 


D5a 


12.0 


2.3 


0.040 




0.79 


D5b 


12.0 


7.1 


5.0 X 10-3 




0.13 


D5c 


12.0 


14.1 


4.0 X 10-3 




0.025 


E14b 


27.5 


9.4 


0.010 




0.010 



Table 3. The mass corresponding to the peak of the mass dis- 
tributions measured in Fig. 1101 and the corresponding mass pre- 
dicted from the Padoan-Nordlund model using our initial condi- 
tions. 



There are two further limitations to our calculation of 
the star formation efficiency that arise from the numeri- 
cal limitations of our model. First of all, many of the cores 
will continue to accrete gas after the time the simulations 
stopped due to violation of the Truelove criterion. This will 
lead to an increase in the estimated values of the SFE. Fur- 
thermore, our periodic boundary conditions artificially pre- 
vent scattering of cores out of the cl o ud an d enhance the 
likelihood of interactions (e.g. lLi et al.l (|2004l )). Allowing the 
cores to disperse away from one another into lower-density 
g as would serve to further decrease the SFE (e.g. the work 
of iBate et all l|2002l )). 



5 DISTRIBUTION FUNCTIONS OF 
FLUCTUATIONS AND CORES 

In this section, we provide the complete set of statistical 
properties of cores that form in our simulations. These in- 
clude the distribution of core masses, angular momenta, 
spins (fl), magnetizations (/3), radii, and mean densities. 
These distribution functions are calculated at the end of our 
simulations, which are determined by the times at which we 
can no longer adequately follow the collapse of the densest 
core due to the violation of the Truelove criterion (ie, the lo- 
cal Jeans length can no longer be adequately resolved). The 
precise time at which this occurs varies from simulation to 
simulation. 



5.1 Mass Distribution 

We have plotted the distribution of masses for all of the 
condensations within our simulations in Fig. (I1U|I . In all 
cases, we have plotted the set of all condensations with a 
dashed line, and the gravitationally bound cores only in 
solid lines. The bound core distribution generally consists 
of a power-law behaviour at high masses, with a turnover at 
10"'' — 10~^mj (when scaled to the initial values in Table 
[T] this corresponds to 0.1 — OArriQ). This type of behaviour 
is c onsistent with other results reported in the l iterature 
Gammie et all (|2003l ) has a turnover at O-Sm^; 



(20041' 



Li et al 



mj 



; iBate fc Bonnell 



find a turnover at ~ 1 — 2a; 10 
20051 ) has the turnover at 0.07mj), and consistent with ob- 



iLuhman et aLll2000l . |2003| ). We produce bound cores with 
masses an order of magnitude below this peak, and thus be- 
low the hydrogen-burning limit,^as suggested should be the 
case for supersonic turbulence (jPadoan fc Nordlund 2004, ). 

We attempted to fit the bound core mass distribu- 
tion w ith the theoretical distribution of lPadoan fc Nordlundl 
This distribution posits that the high-mass part of 
the spectrum is a power law function that arises from the ori- 
gin of the cores in shocks. There is a turnover at low masses 
due to the fact that even the shocks cannot compress small 
enough fluctuations to the point where they become Jeans 
unstable. The Padoan-Nordlund distribution has the form 



dN{m) 
dlog m 



erf 



V2 



ln(m) + %/2o 



(14) 



servations of low-mass star forming regions l|Luhmanll2OO0l : 



where the power-law -1.29 exponent is derived from the ex- 
ponent 7 of the turbulent power spectr um (in our case, 7 = 
—5/3) , via the equation —3/(4-1-7) (see lPadoan fc Nordlundl 
(j2002l )). This exponent arises from the density jump in a 
magnetized, isothermal shock, and an assumption that the 
number of cores in a given volume scales in a self-similar 
manner with the volume, ct is a measure of the width of the 
density probab ility distribution function; i n the turbulence 
model used bv IPadoan fc Nordlundl (|2002[ ). it is related to 
the Alfven Mach number via 0-2 ~ ln(l + \M\). The er- 
ror function arises from the requirement that the mass of 
any given turbulence fluctuation be greater than the local 
Jeans mass in order for collapse to proceed, and leads to a 
turnover at low masses. As the likelihood of a fluctuation 
being sufficiently compressed in a shock so that it becomes 
Jeans unstable is greater if the turbulence is stronger (and 
thus the stronger shocks can compress the fluid more), the 
width of this distribution is directly related to the Mach 
number of the turbulence. We flt (dotted line in Fig. I10|l 
this distribution to our bound core mass function, allowing 
the local Mach number and the height to vary. For compar- 
ison, we also provide the predicted spectrum of the PN02 
model using the initial Alfven Mach number - a curve that 
is indicated by a dot-dashed line. 

We flnd that in most of our simulations, the shape of 
the distribution of bound core masses is reproduced well by 
the Padoan-Nordlund spectrum. As is readily seen by com- 
paring these curves, their model does not fit the peak mass 
of the distributions with one exception (E14b), as demon- 
strated in Table [3] The greatest discrepancy is that the po- 
sition of the peak of the mass spectrum in our simulations 
is independent of the initial Alfven Mach number, contrary 
to the predictions of Padoan & Nordlund. We also find that 
the distribution of cores produces far more low-mass cores 
than would be predicted by the Padoan-Nordlund model for 
our initial Alfven Mach numbers. The trend towards lower 
masses could in part be explained by the important role 
played by the surface pressure in binding the cores (see Fig. 
[8|, which could allow smaller fluctuations to remain over- 
dense long enough for gravitational collapse to begin. 

In the non-magnetic simulations of TP04, as soon as 
one core began to collapse, it collapsed to the point that the 
simulation could no longer follow it due to violation of the 
Jeans criterion. It did so before any other core was able to 
evolve to a signiflcant degree, and thus there was always one 
high-mass outlier in the mass distributions. The same effect 
can be seen in these simulations that are highly supercriti- 
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Figure 10. Distribution of core masses at the end of each of the simulations. The distribution marked by a soUd line contains only 
those cores that are gravitationally bound; the distribution of all cores is marked by a dashed line. The dot-dashed line is the predicted 
spectrum of Padoan & Nordlund, using our initial conditions; the dotted line is a fit to the same spectrum to the bound core distribution. 



cal; the first core to collapse does so quickly enough that it 
alone accretes a significant portion of the mass before the 
simulation is halted. In the runs with a stronger magnetic 
field (and are thus less supercritical), this effect is much less 
pronounced and can even disappear completely. 

Our initial conditions are similar to that of iLi et al.l 

l|2004 ). who find a similar result for the mass spectrum of 
bound cores. As these authors utilize a more restrictive def- 
inition of a core (in that they don't allow for fluctuations 
that are bound by virtue of their surface pressure), the cor- 



respondence between our results and theirs suggests that we 
are not biasing o ur res ults by using a broader definition of 
a core. iLi et al.l (|2004l ') examine the evolution of the mass 
distribution and find that it becomes shallower with time 
as cores merge. We expect to find a similar pattern here, 
although we do not examine it further. 
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Figure 10 - continued 



5.2 Angular Momentum Distribution 

Fig. plots a histogram of the distribution of specific an- 
gular momenta of the cores in each simulation. As in Fig. 
HlOp . the solid lines mark the distribution for the bound 
cores, and the dashed lines the distribution of all fluctua- 
tions, bound and unbound. There does not appear to be 
any clear relationship between the peak of the angular mo- 
mentum distributions of the cores in Fig. with either 
nj or /3. In all of the simulations, the core angular mo- 
mentum distribution peaks at ~ 5 — lOcsL, while the un- 
bound condensation angular momentum distribution peaks 
at the slightly lower values of ~ 1 — ScsL. (For compari- 
son, CaL ^ W^'^ — 10^^ cm^s~^ for our initial conditions.) 
T his is comparab l e to the results of ~ 4 x 10^^ cm^s~^ 
of iGammie et all (|2003|). and significantly larg e r tha n the 
2 X 10^^ cm^s"^ that iBurkert fc Bodenheimerl (|2000h find 
from simulated observations of turbulent noise in an oth- 
erwise static core. It is also several orders of magnitude 
la rger than the angul a r mom enta of prestellar cores found 
in ljappsen fc KlessenI (|2004l 'l. although in that case the dif- 
ference is likely due to the fact that they define cores as col- 
lapsed objects represented by sink cells, and thus exist on 
a much smaller scale than what we can resolve. Observed 
velocity gradients in molecular cloud cores have been inter- 
preted as rotation, with estimate specific angular momenta 
of i ^ 10^^-10^^ cm"^ ijGoodman et al.ll 19931 : ICaselU et al.1 
l2002al ). slightly less than our estimates. We note that the 



Run 


/3 of Simulation 


Median /3 of cores 


B5b 


0.9 


0.34 


B5c 


5.0 


3.2 


B5d 


10.0 


10.9 


C5c 


5.0 


5.3 


C5d 


10.0 


11.4 


C5e 


50.0 


53.9 


D5a 


0.1 


0.090 


D5b 


1.0 


0.87 


D5c 


3.9 


4.5 


E14b 


1.0 


0.91 



Table 4. Initial /3 of each simulation and the median /3 of the 
cores produced in each simulation. 



specific angular momenta of our cores decreases with radius, 
which we will describe in more detail in a future paper. 

We see a bimodal angular momentum distribution for 
run B5b, similar to the bimodal mass distribution. There 
is not a clear explanation for this, although it too is likely 
related to the differences in the mechanics of the collapse in 
a sheet versus collapse in the envelope. 
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Figure 11. Distribution of core specific angular momenta. 



5.3 Magnetization (/3) and Mass-to-Flux (F) 
Distributions 

The mean magnetization ( /3) of each core for our standard 
set of simulations is plotted in Fig. (|12|) . As in Fig. (|10|) . 
the solid line marks the distribution for the bound cores 
that were identified in that simulation, and the dashed line 
marks the distribution of all fluctuations. The median value 
of l3 for the bound cores is within 15% of the original value 
of /3 in all of the simulations except B5b and B5c (see Table 
[4|; in B5b we see a large number of cores with very small 
values of /J that is separate from the distribution centred 
on the mean /9 of 0.9, and similarly in B5c we see that the 
distribution of /3 is strongly skewed towards small values, 
unhke what we see in the other simulations. 



The distribution of the mass-to-flux ratio of the cores 
is plotted in Fig. (|13p . with solid lines denoting bound 
cores and dashed lines denoting all fluctuations. We see that 
nearly all of the bound cores that are produced have mass-to- 
flux ratios less than the initial mass-to-flux ratio of the origi- 
nal simulation domain. We also find that significant numbers 
of subcritical cores can be produced, even in the simulations 
that start in a highly supercritical state. This important re- 
sult s uggests that the observed near-critical mass-to-flux ra- 
tios of lCrutcheij l| 1999 1) can develop from a cloud that is ini- 
tially supercritical. Thus, it is possible to avoid the sheet-like 
collapse behaviour of an initially near-critical cloud while 
still producing near-critical cores similar to those seen in 
star- forming regions. 

In the low-P runs, the separation between the cores in 
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Figure 11 - continued 



the high-density sheet and low-density atmosphere is espe- 
cially pronounced in Figs. (|12I13|I . As the collapse into these 
sheets preferentially occurs along field lines, the magnetic 
field strength both inside and outside of the sheet is ap- 
proximately constant; since /3 oc p/ , the density contrast 
between sheet and environment is directly reflected in the 
/3 distribution for these marginally supercritical runs. It is 
apparent that in these marginally supercritical simulations, 
the magnetic field is not significantly compressed. 

5.4 iltff Distribution — Upper Limits on the 
Binary Frequency? 

A quantity that has proved useful in determining whether 
a collapsing core fragments is the product of O, the mean 
angular rotation rate of the core, and fg, the free-fall 
time at the centre of the core (|Matsumoto fc Hanawair2003l : 
iBaneriee et al.|[2004l ). We calculate this product using the 
specific angular momentum, central density and mean ra- 
dius of each condensation: 

The results are plotted in Fig. (O. 

For most of the simulations, Q.tfi is around 0.1. 
The collapse of such object s will produce discs , which 
iMatsumoto fc Hanawal (|2003l 'l and IBaneriee et~al] (|2004l ') 
suggest will fragment into a combination of rings and bars, 
depending on the amount of anisotropy. Cores with Q,tfi < 



0.03 tend to collapse smoothly into discs without significant 
fragmentation, but very few of our cores fall into this cat- 
egory. Cores with Qtg > 0.3 will not collapse beyond some 
minimum radius, where rotation will cause the core to ex- 
pand again. About 20% of the cores in our simulations fall 
into this last category. We see some notable exceptions for 
simulations with low P. For Run B5b, fltg tends to be higher 
as the sheet that forms produces a number of cores with very 
low central densities and masses away from the sheet; the net 
result is a significantly larger free-fall time. Conversely, run 
D5a has extremely low values of i}tft as the bound cores gen- 
erally have extremely large radii (Fig. llSp . thus resulting in 
a low mean rotation rate. These distributions of Qta there- 
fore say something very interesting about the binary frac- 
tion expected in young star clusters. Our d ata suggest that 
it sho uld b e rather high, in agreement w ith lBate fc BonneUl 
l|2005l ) and lDelgado-Donate et al] l|2004) . 



5.5 Radius and Central Density Distributions 

In Fig. [15] we plot the distribution of the radii of the cores 
found in our simulations, scaled to the default values given 
in Table [T] In general, the bound core radius distribution 
peaks at 10% of the box size of the simulation. As a result, 
most of the cores produced by these simulations are 0.02-0.2 
pc in s ize. This is t h e sam e scale as where most of the cores 
in the Ijiiina et al.l (|l999l ) dataset are found. Only a very 
tiny fraction of cores, ~ 1%, are resolved with fewer than 
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Figure 12. Distribution of mean core /3. 



10 zones. In contrast, 10-20% are typically resolved with at 
least 30 zones. 

The distribution of the peak density of each core is pre- 
sented in Fig. 1161 The existence of a core that is undergoing 
runaway accretion is evident in all of the highly supercriti- 
cal simulations, as this core has a significantly higher den- 
sity than any of the rest. This is less of an effect in the 
stronger-field simulations, due to the reduced compressibil- 
ity of the gas from the presence of the strong magnetic fields 
llPassot et al.[ ,1995: Balsara et al . 2001; H eitsch ct al. 200?; 
IVazguez-Semadeni et alj|2005bi ). The peak of these distri- 
butions are between 10* - lO'^ cm , densities that are well- 
traced by ammonia maps (although it should be noted, we 
chose our initial conditions so that this would be true) and 



again agree with ljiiina et al.l l)l999l l who characterize bound 
cores through ammonia surveys. 



6 DISCUSSION AND CONCLUSIONS 

We have presented a detailed examination of the physical 
properties of cluster forming cores as they arise in sim- 
ulations of decaying turbulence in magnetized, turbulent 
"clumps" within molecular clouds. This effort is aimed at 
trying to understand the origin of the core mass function, 
as well as many other properties of cores and eventually, the 
stars that form within them. In pursuing this goal, we have 
been careful to examine all of the physical forces that are 
acting upon the fluctuations in our simulations - our virial 
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analysis of cores is the most exhaustive analysis we can find 
in the literature. Our results are also not sensitive to the na- 
ture of the initial spectrum of the velocity field fluctuations 
that is imposed on the simulations. 

We find that an initial turbulent, magnetized molecu- 
lar cloud must be significantly supercritical if its subsequent 
fragmentation is to resemble at all the properties of evolved 
molecular cloud cores that we see today. The cores that are 
produced in these highly supercritical simulations can still 
have (local) mass-to-fiux ratio s that are close to critical, re- 
producing the trends seen bv lCrutcherl l|l999D for observed 
cores. Simulations that are not initially significantly super- 
critical preferentially collapse along field lines to form sheets, 
and the fragmentation of these sheets produce much flatter 
mass distributions. 

We flnd through the use of the virial theorem 
that the surface pressure and surface magnetic field 
play a critical role in creati ng bound cores (see also 
iBallesteros-Paredes et ai] ([l99i)). 

We obtain a distribution of masses for the bound cores 
that shows a power-law at large masses for simulations that 
began with moderately supercritical magnetic fields. The 
peak of these mass distributions are robustly found at 0.001- 
0.01 times the initial Jeans mass of the simulation. We fit 
these cor e mass distributions with th e theoretical mass spec- 
trum of |Padoan"fc^ Nordlundl (|2002l ) and generally obtain 
reasonable fits to the shape of the core mass spectra us- 
ing their model (albeit with fewer very-low-mass cores pro- 
duced than expected by the Padoan-Nordlund mass spec- 



trum). Cores that are only marginally supercritical have a 
significantly flatter mass distribution, and do not appear to 
develop a clear power-law at large masses. However, the peak 
mass of the core mass spectra in our simulations are not well 
flt by this model. In particular, there is little sensitivity of 
the peak mass in our simulations to the initial Mach number 
of the turbulence as proposed by their model. 

The bound cores that are produced have specific an- 
gular momenta that, when scaled to appropriate physical 
units, are on the order of 10^^ — 10'^'^ cm^s~^. On the ba- 
sis o£_thsse_S2m_£ates_aBd the d etail ed collapse calculation s 
of lMatsumoto fc Hanawal (|2003h and lBaneriee et al.l (|2004h . 
many of these cores are likely to further fragment during 
the course of their collapse, likely leading to the formation 
of binary or multiple systems. 

We find that the bound cores have a distribution of mag- 
netization (/3) that has a median value that is approximately 
the same as the mean f3 of the simulation. The mass-to-flux 
ratio (r) of these cores, however, is generally less than the 
mass-to-flux ratio of the original simulation, by up to an or- 
der of magnitude. This suggests that a distribution of cores 
can be produced that have mass-to-flu x ratios that are close 
to critical (as seen by ICrutcheilll 19991 ) ) from a cloud that is 
initially highly supercritical. 

The distribution of core radii is a narrowly peaked func- 
tion centered at 0.02-0.2 pc. The distribution of core den- 
sities peaks at the mean density of the simulation, but is 
skewed towards higher densities due to gravitational col- 
lapse. The marginally supercritical runs develop a double 
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Figure 13. Distribution of the mass-to-flux ratio T of cores. 



peak in their density distributions, with the one peak aris- 
ing from the cores collapsing in a dense sheet and the other 
from cores in the low density surroundings. 

Wc conclude that simulations of magnetized, decaying 
turbulence can provide an excellent account of the core mass 
functions and a wide range of core properties in such clouds. 
The role of magnetic fields can indeed be very significant in 
the dynamics of molecular clouds and their cluster forming 
cores. Indeed our best fits to a wide variety of observations of 
cloud cores suggest that the strong magnetization of cores 
that is often observed is the result of compression by the 
turbulence, and is not characteristic of the entire volume of 



the cloud. This further suggests that turbulence, and not 
wide spread cloud magnetic field, lies at the heart of the 
origin of core formation and the origin of the IMF. 
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Figure 16. Distribution of the peak core density, scaled to the initial values in Table [Tl 
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